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Abstract 

We report a neutron scattering study of the spin-1/2 alternating bond anti- 
ferromagnet Cu(N03)2 • 2.5D2O for 0.06 < ksT/Ji < 1.5. For ksT/Ji < 1 
the excitation spectrum is dominated by a coherent singlet-triplet mode cen- 
tered at Ji = 0.442(2) meV with sinusoidal dispersion and a bandwidth of 
J2 = 0.106(2) meV. A complete description of the zero temperature contribu- 
tion to the scattering function from this mode is provided by the Single Mode 
Approximation. At finite temperatures we observe exponentially activated 
band narrowing and damping. The relaxation rate is thermally activated and 

wave vector dependent with the periodicity of the reciprocal lattice. 
75.10.Jm, 75.40.Gb, 75.50.Ee 
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Transverse phonons and spin waves are propagating small amplitude oscillations of a 
static order parameter in a broken symmetry phase. Isotropic quantum antiferromagnets 
with a gap in their excitation spectra can also support coherent wave-like excitations, but 
these differ from phonons and spin waves in that they move through a system with no static 
order. Specific examples of such systems include the spin-1 chain []T|, even- leg ladders |0] 
and the alternating bond spin- 1/2 chain Since they are not based on the existence 

of a static order parameter that sets in at a well defined transition temperature, coherent 
excitations in these systems are expected to emerge smoothly with decreasing temperature 
as short-range correlations develop. 

In this letter we document this unique cooperative behavior through an experimen- 
tal study of the temperature dependence of magnetic excitations in an isotropic, order- 
parameter- free quantum magnet. Specifically, we have studied magnetic excitations in the 
alternating spin-1/2 chain Copper Nitrate (CN) as a function of wave vector, energy, and 
temperature. The spin Hamiltonian for this system can be written 

2n+l ■ ^2n+2 )• (1) 

n 

Because J2/ Ji ~ 0.24 is small it is useful to think of CN as a chain of pairs of spins- 

1/2. Each pair has a singlet ground state separated from a triplet at Ji ~ 0.44 meV. The 
weak inter-dimer coupling (J2 ~ 0.11 meV) yields a collective singlet ground state at low 
temperatures with triplet excitations that propagate coherently along the chain. We have 
characterized dynamic spin correlations in the temperature range 0.06 < ksT / Ji < 1.5, in 
considerable detail. We find that heating yields thermally activated band narrowing and an 
increased relaxation rate that varies with wave vector transfer with the periodicity of the 
reciprocal lattice. 

CN (Cu(N03)2 ■ 2.5D2O), is monoclinic § (space group I12/cl) with low T lattice 
parameters a = 16.1 A , 6 = 4.9 A , c = 15.8 A , and (3 = 92.9°. The vector connecting 
dimers center to center is uq = [lll]/2 for half the chains, and Uq = [lIl]/2 for the other 
half. The corresponding intra-dimer vectors are di = [0.252, ±0.027, 0.228] respectively. 
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In our experiment, the wave vector transfer Q was perpendicular to b so the two sets of 
chains contributed equally to magnetic neutron scattering. The sample consisted of four 
92% deuterated, co-aligned single crystals with total mass 14.09 g [Q. 

Inelastic neutron scattering measurements were performed on the inverse geometry time 
of flight spectrometer IRIS at the Rutherford Appleton Laboratory, UK [^. Disk choppers 
selected an incident spectrum from 1.65 meV to 3.25 meV pulsed at 50 Hz, and a backscat- 
tering pyrolytic graphite analyzer bank selected a final energy Ef = 1.847 meV. The Half 
Width at Half Maximum (HWHM) elastic energy resolution was 10.5 /ieV. The b direction 
of CN was perpendicular to the horizontal scattering plane and [101] pointed towards the 
low angle part of the analyzer bank at an angle of = 20(1)° to the direct beam. The 
direction for rotation of the a axis into the c axis coincided with the direction of decreasing 
scattering angle. In this configuration the projection of wave vector transfer on the chain 
Q\\ = ki cos (p — kfCos{(f) — 26) takes on a unique value for each detector in the range of 
scattering angles 20° < 26 < 160° covered. We can therefore present our data as a function 
of energy transfer hu and wave vector transfer along the chain g = Q ■ Uq. In addition there 
is a specific value of Q± = ki sin0 — kf sin(0 — 26) associated with each (g, hu) point such 
that sensitivity to dispersion perpendicular to the chain is maintained in the projection. 
Count-rates were normalized to incoherent elastic scattering from the sample to provide 



absolute measurements of I{Q,uj) = \^F{Q)\'^2S{Q,iu). Here g = ^ {gl + g\)/2 = 2.22 [| 
F[Q) is the magnetic form factor for Cu^"*" [|ll|, and S{Cl,uj) is the scattering function [jl2 



Figures |l|(a)-(c) show the measured neutron scattering spectrum at T = 0.3 K, 2 K, and 4 
K. Focusing at first on the 0.3 K data, we observe a resonant mode centered around Ji = 0.44 
meV with bandwidth ^ J2, consistent with the predictions of perturbation theory [|13|. The 



mode energy has the periodicity 27r of the one-dimensional reciprocal lattice with minima 
for q = Q ■ uq = n ■ 2tt, indicating antiferromagnetic inter-dimer interactions. The intensity 
of the mode varies with a periodicity that is incommensurate with that of the dispersion 
relation. There is a simple explanation for this, namely that the intra-dimer spacing that 
enters in the neutron scattering cross section is incommensurate with the period of the 



alternating spin chain [|^. An exact sum-rule for iS(Q,C(j) provides the following direct link 



between the microscopic structure and the intensity of inelastic neutron scattering |T5 

/oo 
ujS{Q„ uj)duj 
-oo 



= -^^Jd(So-Sd)(l-cosQ-d), (2) 

where {d} is the set of all bond vectors connecting a spin to its neighbors. Figure 2 (a) 
shows h{uj)Q at T = 0.3 K and T = 4 K derived by integrating the corresponding data sets 
in Fig. 1. The solid lines show fits based on Eq. |^ including only intra-dimer correlations 
(So ■ Sdi) and a constant to account for multiple scattering (see below). The excellent 
agreement between model and data provides direct evidence for singlet formation between 
spins separated by di. The modulation amplitude is Ji(So ■ Sdi). 

From the wave vector dependence of energy integrated intensities we turn to spectra at 
fixed wave vector transfer. Figure 3 shows cuts through raw data for q = 2tt and q = Sir 
at three temperatures. At T = 0.3 K we see resolution limited peaks. Upon increasing 
temperature these peaks broaden and shift towards the center of the band. Careful inspection 
also reveals that the higher energy peak at g = 37r broadens less than the lower energy peak 
at g = 2n. From gaussian fits to constant energy cuts such as these, we extract the resonance 
energy and half width at half maximum versus wave vector transfer shown in Fig. 2 (c) and 
(d). This analysis shows that the finite T relaxation rate is wave vector dependent with the 
apparent periodicity of the reciprocal lattice. 

We now proceed to extract the detailed temperature dependence of the integrated in- 
tensity, the effective bandwith and the relaxation rate. To take full advantage of the wide 
sampling of Q — space in our experiment and to account for resolution effects, the analysis 
is based on "global fits" of the following phenomenological form for iS(Q, u) to the complete 
Q and u dependent data set at each temperature. For hu > we write 

5(Q, c.) = h^JT^Ji^^ - <Q))> (3) 

e(Q) 1 - exp(-/?e(Q)) 

where f{E) is a normalized spectral function. The other terms implement the first moment 



sum-rule of Eq. Q in the limit where f{E) is sharply peaked on the scale of Ji. Eq. |^ 
represents the "Single Mode Approximation" (SMA) that has been used with success to link 
the equal time structure factor and the dispersion relation for collective modes in numerous 



many body systems [|T^-|T^. Given that two magnon scattering carries less than 1% of the 



spectral weight for CN |T^, the SMA should be excellent at sufficiently low T. For the 
dispersion relation we use the following variational form based on first order perturbation 
theory H 

e(Q) = Ji - JuCosQ ■ u, (4) 

^ u 

The vectors {u} connect neighboring dimers center to center, both within and between the 
chains. For comparing to the experimental data, iS(Q, u) was convolved with the instrumen- 



tal resolution [|T0|. In all fits it was necessary to take into account multiple neutron scattering 
events involving elastic incoherent nuclear scattering followed or preceded by coherent in- 
elastic magnetic scattering. We believe that such processes are responsible for the weak 
horizontal band of intensity that is visible in Fig. |l](a). For T = 0.3 K we used f{E) = 6{E) 
and obtained the fit shown in Fig 1 (d). To better evaluate the quality of the global fit we 
also show cuts through the model calculation as solid lines in Fig. 3 (a). There is excellent 
agreement between model and data with an overall prefactor and four exchange constants 
as the only fit parameters. The prefactor refined to (Sq • Sdi) = —0.9(2) consistent with 
the value of -3/4 expected for isolated singlets. Allowing for inter-dimer correlations by 
including the corresponding term from Eq. || in the global fit yields (So ■ Sda) = —0.04(8). 
From the best fit parameters in the variational dispersion relation we get Ji = 0.442(2) meV 
and J2 = 0.106(2) meV for the two intra-chain interactions and Jl = 0.012(2) meV and 
Jr = 0.018(2) meV for dimers separated by = [|,0,0] and = [0,0, |] respectively. 
Because b ± Q throughout the experiment our data do not yield an estimate for inter-chain 
interactions between dimers displaced by [1 =F ll]/2. However, in a separate measurement 



with the crystal oriented in the (hkh) plane we were able to place an upper limit of 0.02 
meV on the corresponding parameter in Eq. ^. The value J2/J1 = 0.240(5) that we obtain 



is measurably smaller than the value J2/J1 ~ 0.27 derived from magnetic susceptibility 
data. A likely explanation for this discrepancy is that susceptibility measurements cannot 
distinguish intra-chain from inter-chain interactions. 

To analyze the finite T data we replaced the spectral function with a normalized Gaussian 
with HWHM 

Y 

r(g) = To + y cosg, (5) 

The functional form for the dispersion relation (Eq. ^) was maintained, but to account 
for the bandwidth narrowing that is apparent in the raw data, we introduced an overall 
renormalization parameter relating finite T "effective" exchange constants in the variational 
dispersion relation to the bare temperature independent exchange parameters: = n(T)J■^^. 
The fits obtained are excellent as can be ascertained by comparing the left and right columns 
in Fig. |l] and the solid hues through the data in Fig. |^. The sohd lines in Fig. ^(b) and (c) 
show the dispersion relation and g-dependent relaxation rate derived from the global fits. 
They are consistent with the data points derived from constant-g cuts indicating that the 
variational forms employed for e(g) and T{q) have not biased the global fitting analysis. We 
note that the quality of the fits did not change significantly when using a Lorentzian rather 
than a Gaussian spectral function. So while the data provide reliable information on the q- 
and T-dependent relaxation rate, they do not accurately determine the spectral function. 

The temperature-dependent parameters derived from this analysis are shown in Fig. ^. 
The prefactor for global fits at each temperature yields the intra-dimer spin correlation 
function, which we plot versus T in Fig. §(a). As expected |(So • Sdi)| decreases with 
increasing T as the populations of the four states of each spin pair equalize. For an isolated 
spin pair it can be shown that (So-Sd) = —(3/4)An(/3Ji), where An{j3Ji) = (1 — e~'^'^i)/(l + 
3g-/3Ji^ is the singlet triplet population difference. After fitting a scale factor, this form 
provides an excellent description of the temperature dependence of the data (solid line in 
Fig. ^(a)). The Random Phase Approximation (RPA) applied to interacting spin dimers 
p0|^3| predicts that the bandwidth renormalization factor, n{T), also follows the singlet 
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triplet population difference. A similar result holds when the singlet ground state is induced 
by single ion anisotropy |2^. We compare n{T) to An{i3Ji) in Fig. ^(b). While there is 
qualitative agreement, the RPA clearly predicts more bandwidth narrowing than is actually 
observed. 

Triplet relaxation is due to scattering from the thermal ensemble of excited states. Be- 
cause the density of triplets is thermally activated, the relaxation rates should be too. Fits 
of a simple activated form Fj = 7jexp(— Aj/ZcfiT) to the T < 4 K data in Figs. §(c) and (d) 
give Ao = 0.24(2), 70 = 0.10(2) meV, Ai = 0.32(3), and 71 = 0.08(2) meV. If we instead 
fix Aj = Ji and allow for a power law prefactor : Fj(T) = 7j( Ji/fc^T)"' exp(— Ji/Zcb?") we 
also obtain excellent fits with = 1.0(2), 70 = 0.13(4) meV, ai = 0.6(2), and 71 = 0.10(1) 
meV. 

In summary, we have examined dynamic spin correlations in the strongly alternating 
spin chain Cu(N03)2 ■ 2.5D2O for 0.06 < ksT / Ji < 1.5. For fc^T ^ Ji we find a coherent 
dispersive triplet mode whose contribution to the scattering function is perfectly accounted 
for by the SMA, and we have determined accurate values for inter-dimer interactions in the 
material. Upon heating, the sharp dispersive mode gradually deteriorates through band- 
narrowing and the development of a wave-vector dependent lifetime. A semiclassical theory 
for finite-T excitations in gapped spin chains was recently developed by Sachdev and Damle 



| 25| . It relies on A/ J being a small parameter as is the case for the Haldane phase of spin-1 
chains, and in weakly dimerized spin-1/2 chains. In CN the spin gap is instead much greater 
than the magnetic bandwidth. Our data should provide a focus for theoretical attempts to 
describe finite temperature properties in this "strong coupling" limit. 
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FIGURES 

FIG. 1. Normalized scattering intensity, I{q, u;), for Cu(N03)2 -2.5020 at T = 0.3 K (a), T = 2 
K (b), and T = 4 K (c). Maximum intensity on the color scale shown above are Imax = 28 meV~^, 
16 meV~^ and 8 meV^^ respectively, (d)-(f) show model calculations based on Eq. ^ 

FIG. 2. First moment 

computed from Fig. 1 (a) and (c) using h{uj)Q = fi^ uj{1 — e~^^'^)S{Q,uj)duj (a), resonance 
energy (b), and HWHM SE{q) (c) of the triplet mode in copper nitrate at T = 0.3 K and T = 4 
K. The points in (b) and (c) were obtained from fits to constant- q cuts through the data in Fig. 
1. Dashed line in (c) shows the HWHM instrumental energy resolution at hoj = 0.45 meV. Solid 
lines in (b) and (c) were obtained from global fits at each T. 

FIG. 3. Constant-g cuts through the data shown in Fig. 1 at g = 27r and Stt. The solid lines 
are the results of global fits at each T. Dashed lines indicate peak positions at T = 0.3 K, and 
highlight band narrowing at higher T. 

FIG. 4. Parameters characterizing the temperature dependence of dynamic spin correlations in 
CN, obtained from fits shown in Figs. 1 and 3. The solid lines in (a) and (b) show RPA theory. 
The solid lines in (c) and (d) are fits to Ti{T) = 7j( Ji/Zc^T)"* exp(— Ji/Zc^T). 
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http://arXiv.org/ps/cond-mat/9908201v2 
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